CrossMark 



research papers 



Journal of 

Applied 

Crystallography 

ISSN 0021-8898 



Received 18 October 2012 
Accepted 12 December 2012 



Directional pair distribution function for diffraction 
line profile analysis of atomistic models 

Alberto Leonard!^ Matteo Leoni and Paolo Scardi* 

Department of Civil, Environmental and Mechanical Engineering, University of Trento, Via Mesiano 
77, Trento (TN), 1-38123, Italy. Correspondence e-mail: paolo.scardi@ing.unitn.it 



The concept of the directional pair distribution function is proposed to describe 
Hne broadening effects in powder patterns calculated from atomistic models of 
nano-polycrystalline microstructures. The approach provides at the same time a 
description of the size effect for domains of any shape and a detailed 
explanation of the strain effect caused by the local atomic displacement. The 
latter is discussed in terms of different strain types, also accounting for strain 
field anisotropy and grain boundary effects. The results can in addition be 
directly read in terms of traditional line profile analysis, such as that based on the 
Warren- Aver bach method. 



1. Introduction 

Considerable efforts in the study of nanomaterials focus on 
the visualization of strain at the atomic level. Nevertheless, 
even the most advanced experimental methods based on 
coherent diffraction are confined to single isolated nano- 
particles and relatively weak strain fields (Thomas, 2008; 
Robinson & Harder, 2009; Watari et a/., 2011). Strains in nano- 
polycrystalline microstructures are still a challenge to direct 
visualization techniques, because of the variety of strain 
sources and the level and complexity of their distribution. 
In this context, diffraction line profile analysis (LPA), 
although indirect and much less detailed than coherent 
diffraction and other microscopies, is still one of the most 
convenient and widely used experimental techniques. Further 
insights into the mechanisms and distributions of strains can 
be produced by atomistic models, in particular by molecular 
dynamics (MD) simulations (Van Swygenhoven et al, 2000; 
Van Swygenhoven, 2002; Derlet et aL, 2005; Derlet & Van 
Swygenhoven, 2006; Cao et al, 2008; Stukowski et al, 2009). So 
far, however, the connection with diffraction LPA has been 
discussed in terms of simplified integral breadth methods only 
(Derlet et al, 2005; Leonardi et al, 2011): the potential of 
joining MD simulations with state of the art LPA is still largely 
unexplored. 

In the past decade LPA evolved in the direction of Fourier 
methods based on physical models of the microstructure 
(Mittemeijer & Scardi, 2004). Despite success in many appli- 
cations, even the most advanced methods of whole powder 
pattern modelling (Scardi & Leoni, 2002, 2001) still rely on a 
simplified description of lattice defects, and cannot be proven 
to be univocal in the identification and quantification of the 
different effects contributing to the line profiles. For example, 
the uniqueness of the anisotropic line broadening due to the 
strain field of dislocations, the effects of grain-to-grain defor- 
mation fields and the role of grain boundaries are currently 
open questions. It is not clear to what extent these different 



sources of lattice deformation can be identified by even the 
most advanced LPA methods. 

MD simulations of nano-polycrystalline microstructures can 
give useful insights into these problems, provided that the 
usual concept of strain is reconsidered and related in a more 
direct way to diffraction LPA. Inhomogeneous strains have a 
pecuHar effect on diffraction, which is sensitive to the distri- 
bution of the strain component projected along the scattering 
vector direction; this strain component needs then to be 
considered as a function of the separation distance between all 
possible pairs of atoms (or unit cells, on a coarser size scale) in 
each given crystalline domain (Warren, 1990). 

As shown by several MD studies, the Debye scattering 
equation (DSE) is the most direct way to simulate the powder 
diffraction pattern from a given atomistic model of nano- 
polycrystalline microstructure (Cervellino et al, 2003, 2010; 
Derlet et al, 2005; Gelisio et al, 2010, 2011), as it just requires 
atomic positions (strictly speaking, distances between all 
possible pairs of atoms), with no assumption on the crystal 
structure or on the presence of lattice defects. While the DSE 
is rigorous and appealing for simulating the powder pattern, it 
is not so useful for analysing the effects of strain, separating 
different sources of line broadening and studying line aniso- 
tropy, which makes line broadening different for profiles with 
different Miller indices. 

In the present work we introduce the concept of the 
directional pair distribution function (D-PDF) to represent 
the finite size and shape of coherently scattering domains 
(crystallites) and the local atomic displacement in a way 
directly legible in terms of diffraction effects. Starting from an 
atomistic simulation of single crystals as well as of nano- 
polycrystalline microstructures, the D-PDF allows a separa- 
tion of all contributing hkl line profiles in the powder 
diffraction pattern. Each line profile can be described in terms 
of size and strain contributions, in a quite similar way to 
traditional LPA based on Fourier analysis [e.g. the Warren- 
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Averbach (WA) method; Warren & Averbach, 1950, 1952; 
Warren, 1990]. The proposed approach has at least two valu- 
able apphcations: (i) to vahdate different LPA methods and 
better understand their results; (ii) to study MD simulations 
and their relation with real microstructures in terms of a well 
known, easy to perform experimental technique like powder 
diffraction. 

The chosen case of study concerns a metallic nano-poly- 
crystalline system made of randomly oriented grains with 
irregular (though not far from equiaxed) shapes. The system 
was equilibrated (energy minimized and thermalized) by 
conventional MD based on the embedded atom method 
(EAM; Daw & Baskes, 1983), so that no lattice defects other 
than the grain boundaries are present. Under these conditions, 
intergranular strains due to the equilibration process should 
be the only microstructural effects (besides shape and finite 
size of the grains) affecting the line profiles. Detailed infor- 
mation on the direction-dependent strain field and its effect on 
the simulated powder diffraction pattern can be obtained by 
means of the D-PDF concept. The role of the grain boundaries 
is also discussed. 

2. Copper nano-polycrystalline microstructure: 
generation and strain distribution 

A cubic box (side length 260.28 A) was divided into 50 cells by 
the recently developed constrained modified Voronoi tessel- 
lation (CMVT) (Gross & Li, 2002; Xu & Li, 2009; Suzudo & 
Kaburaki, 2009; Leonardi, Scardi & Leoni, 2012; Leonardi et 
al, 2013a). The CMVT allows control of different properties 
during the generation of the cells. In the present case, the 
CMVT parameters were set to produce high sphericity 
(arithmetic mean ^ = 0.800122) and a lognormal distribution 
of grain sizes (for the distribution of diameters of equivalent 
volume spheres: arithmetic mean 82.6 A, standard deviation 
21.4 A); grains were then filled with about 1.5 million Cu 
atoms, using a unit-cell parameter of 3.615 A (Fig. 1). 

Equilibration of the starting ('crystallographic') system was 
done by EAM, using the Cu potential reported by Foiles et al 
(1986). After the energy minimization, an isothermal-isobaric 
time integration at 100 K was performed by the LAMMPS 
software (Plimpton, 1995), reaching equilibrium conditions. 
The latter were assessed by comparing the deformation field 
after 1.2 and 2.4 ns of equilibration times, in terms of both 
volumetric (isotropic) and deviatoric (anisotropic) strain 
components (Fig. 2). The zero strain is referred to an 
equivalent unit cell of [5boxside/(A^atoms/4)]^^^ where ^boxside is 
the equilibrated box side length. 

While the MD simulation was performed on a box with 
periodic boundary conditions (PBCs), grains cut across by the 
cube faces were joined, so as to create a more plausible system 
(Fig. la) for the simulation of the powder diffraction pattern 
by the DSE. 

As shown schematically in Fig. 1(h), strain is not constant 
across the grains: a steep gradient is observed in the outer 
layers, extending from the grain boundary region inward for 
about 10 A. This trend, quahtatively the same for all grains in 



the cluster, is shown in Fig. 2: here the mean volumetric {a) 
and mean deviatoric (b) strain in each grain are plotted 
against the distance from the grain boundary. Insets show that 
the variance of the strain distribution also steeply decreases in 
the sub-surface region. 

While the deviatoric strain decreases to a small constant 
value in the core region, the value reached by the volumetric 
component inside grains depends on the grain size. This 
feature is visible in Fig. 2(a) but it can be seen more directly in 
Fig. 3, where the mean unit-cell parameter of each grain is 
plotted as a function of the diameter of the equivalent volume 
sphere. Larger grains put in tension the smaller ones, with an 
apparently nonlinear trend. In Fig. 3(b) it is also possible to 
see the anisotropy of this deformation, which is higher along 
[/zOO] than along [hhh], respectively corresponding to the 
elastically softer and stiffer directions in copper. 

3. Directional pair distribution function (D-PDF) 

The D-PDF is a convenient concept to represent strain effects 
on the diffraction line profiles. It is obtained by counting the 
number of atom pairs for each distance L along the [hkl] 
direction of the scattering vector in a given grain of the cluster. 
As shown in Fig. 4, the projection of the atom-atom dis- 
placement AL along [hkl] is actually considered: the D-PDF is 




Figure 1 

(a) Copper nano-polycrystalline microstructure after reconstruction of 
grains cut across by the PBCs. (b) Example of volumetric strain along the 
{hOO) directions in a grain. 
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Figure 3 

(a) Type II volumetric strain as a function of the diameter of the equivalent volume sphere; (b) unit-cell parameter as a function of the greatest thickness 
along the direction considered. See text for details. Dot-dashed lines represent the average values. 




Figure 4 

(a) Directional pair distribution function for grain No. 33 in the cluster along the [hOO] direction; (b) example of distribution for L = 14.485 A; the area 
under the curve (shaded) corresponds to the circled value of the normalized common volume function (CVF) in (c). (c) The trend of the normalized area 
of Pi^(AL^qq) as a function of the pair distance. The inset in (c) illustrates the concept of the 'ghost' (see text for details). 
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then made of a series of histogram distributions, centred on 
the mean distance between nth neighbour atoms, ranging from 
the closest neighbour to the longest distance in the grain along 
[hkl]. Each distribution pj^(AL^j^i) can be calculated as 
described above, along [hkl] of a given grain, or be averaged 
over all crystallographically equivalent (hkl) directions in the 
grain, or be averaged over all equivalent directions of all 
grains in the cluster (so-called super D-PDF). 

It is important to note that the D-PDFs along symmetrically 
equivalent directions are in general different, as the strain field 
caused by the neighbouring grains is not subject to any 
symmetry restrictions. These differences, however, tend to 
disappear in the average for a given grain or for the entire 
cluster. 

As an example. Fig. 4(a) shows the D-PDF along the [hOO] 
direction of grain number 33 in the cluster, with one such 
distribution shown in detail in Fig. 4(b). The trend of the 
normalized area of p^(AL^qq) as a function of the pair distance 
in shown in Fig. 4(c). This quantity is equivalent to the 
common volume function (CVF), as introduced by Stokes & 
Wilson (1944), representing the intersection volume between 
a crystalline domain and its 'ghost', i.e. the same domain 
shifted a distance L along [hkl] (inset of Fig. 4c). As shown by 
several authors (Wilson, 1962; Guinier, 1963; Warren, 1990), 
within reasonable approximations the Fourier transform (FT) 
of the CVF gives the so-called size component of the 
diffraction line profile, because of the finite size of the crys- 
talline domain along the given [hkl] direction. As the D-PDF 
is not normalized, the corresponding diffraction line profile 
obtained by FT includes the appropriate weight (volume or 
number of atoms) for the given grain. 

The D-PDF concept allows for a simple and reliable way to 
separate the effects of domain size/shape from those due to 
lattice defects and microstructural features in general. While 
the p^(ALf^,^i) area is related to the domain size/shape 
component of the line profile, the width and shape of the 
D-PDF provide a detailed description of the atomic dis- 
placement over different distances. As shown below, this is the 
same representation of strain effects on line profiles as that 
used by traditional LPA based on Fourier analysis, like the 
WA method (Warren & Averbach, 1950, 1952; Warren, 1955; 
Warren, 1959) and methods proposed by Stokes and Wilson 
(Stokes & Wilson, 1944; Eastabrook & Wilson, 1952). 

It is convenient to introduce the strain over a distance L 
along the [hkl] direction, £hki(L) = AL^j^JL, and the scat- 
tering vector modulus, q — 4n^mO/X, where 0 is half the 
scattering angle and X is the wavelength of the incident 
radiation. As originally shown by Stokes and Wilson (Stokes 
& Wilson, 1944; Wilson, 1962), the FT of the peak profile 
component due to the local atomic strain, y^^(L), can be 
reasonably approximated as 

Y^,,(L) = (exp[/^L^,,,(^)]> = <i(L) + iB^,,(L) 

= f ^w[iqL6hki(L)]pL(^hki) ^^hkh (1) 

where the |Fp term (related to the structure factor, F) is 
omitted for simpHcity. As shown in equation (1), y^/(L) is the 



Fourier transform of the D-PDF, which is a complex quantity 
unless the D-PDF is symmetrical. In equation (1) the D-PDF 
was written for strain, which is straightforward considering 
that PL(£)d6 = pL(AL)L-^dAL. 

The normalization used in equation (1) is such that 
y^^(O) = 1, so it is convenient to represent the CVF in 
normalized form too: 

Al,(L) = N,„{L)/N{Q) = n«(^)/nrain- (2) 

where N^j^i(L) is the number of atom pairs at distance L along 
[hkl] and N(0) is the total number of atoms in a given grain; 
correspondingly, Vf^j^i(L) is the CVF along [hkl] and Vgrain the 
volume of the given grain. It is worth noting here that the 
present description in terms of D-PDF is on a finer and 
generally more precise atomic scale, whereas traditional LPA 
usually considers unit cells as smallest units (Wilson, 1962; 
Warren, 1990). 

Within the limit of Stokes & Wilson's theory (Stokes & 
Wilson, 1944; Eastabrook & Wilson, 1952), the FT of the 
overall line profile for the [hkl] direction in a given grain of the 
cluster is 

= V^,,,,AUL)[A^„(L) + iB^L)]. (3) 

The distribution pj^(AL^f^i) is frequently assumed to be 
symmetrical (5^^-^ = 0) and in particular Gaussian, so that 

Y°,(L) 2. AUL) = exp[-q'L%sl,iL))/2]. (4) 

These are the underlying hypotheses of the WA method 
(Warren & Averbach, 1950; Warren & Averbach, 1952; 
Warren, 1955, 1959, 1990), which can be shown to be stiU 
approximately valid even if the distribution is not exactly 
Gaussian, provided it is a bell-shaped symmetrical distribution 
function (Warren, 1959). As further discussed below, the 
D-PDFs for the studied cluster are neither Gaussian nor 
symmetrical. 

The powder diffraction pattern for a cluster of grains can be 
obtained by adding up all contributions as 

00 

I(q) <x E ^grain E / AI,(L)[AUL) cos(qL) 

grain hkl — oo 

+ Bf„(L)smiqL)]dL, (5) 

where the proportionality factor includes |F|^, several 
constants and known trigonometric terms, whereas the actual 
integration limits are determined by the longest atom-atom 
distance along [hkl] in the given grain. The two sums extend, 
respectively, to all grains and to all crystallographically 
equivalent directions. 

As an example of application of the D-PDF, Figs. 5(a) and 
5(b) show, respectively, the real and imaginary parts of the FT 
for several orders of reflections belonging to the {hOO} family, 
Y^qq(L) = A^qq(L) + iB^QQ(L). The real part decays increas- 
ingly faster for higher orders of reflections, whereas the 
imaginary part sensibly deviates from zero, as an effect of the 
asymmetry of the D-PDFs. The trends of A^qq = AIqqA^qq and 
^hoo = ^hoo^hoo (assuming V^,^^ = 1 for simpHcity) show the 
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combined effect of domain size and strain: the decay of the 
real part is smooth and demonstrates that the integration 
hmits in equation (5) do not need to extend beyond ca 120 A, 
in this specific case, as an effect of the finite size of the domain. 
The fact that different orders of reflections belonging to the 
same family have different trends, with a faster decay for the 
higher orders (i.e. larger q values), clearly demonstrates the 
presence of a strain broadening component. 

4. Powder pattern from a nano-polycrystalline 
microstructure 

As already pointed out in the Introduction, the Debye scat- 
tering equation (Cervellino et al., 2003, 2010; Derlet et al., 
2005) is the most straightforward and correct way to calculate 
the powder diffraction pattern for an atomistic model. The 
DSE makes no assumptions on crystalline structure and lattice 



defects, as it is based only on correlations between afl possible 
pairs of atoms. As such it can be computationally demanding 
(although nowadays entirely viable; Gehsio et al, 2010) but 
quite rigorous, so that we can consider the DSE result as an 
'experimental' pattern, with which the D-PDF analysis 
discussed above can be compared. 

The goal here is to test the D-PDF approach, to better 
understand and vahdate the traditional LPA methods and 
their results, but also to study the features of a nano-poly- 
crystalline microstructure obtained by MD in terms of a well 
known and simple experimental technique like powder 
diffraction. 

As a first test we considered the starting cluster, just after 
the atom filling procedure but before any energy minimization 
and thermalization steps. This corresponds to a system of 
perfect crystalline grains, where the only effect on line profiles 
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q = 47csine/A, (A"') 
Figure 6 

Powder diffraction pattern of the cluster of Fig. 1 as obtained by the DSE (circles) and the D-PDF approach [equation (5), line]. The differences are 
shown below. The inset highlights details in the peak tail region. 
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is that of the finite domain size. As shown in Fig. 6, the match 
between the DSE pattern and that generated by using equa- 
tion (5) is remarkably good (Leonardi et al, 2013c; Leonardi, 
Leoni et al, 2012). This result is to some extent expected, as no 
strain and no deviation from a perfect crystalline order is 
present; however, this step is important to assess the quality of 
the hypotheses underlying equation (5), such as the tangent 
plane approximation (Beyerlein et al, 2011) and the lack of 
grain-grain correlations, which are intrinsically considered by 
the DSE. 

The MD equilibration procedure introduces strains in the 
system, leading to a static component of the atomic displace- 
ment, which adds to the dynamical component due to thermal 
vibrations. At any given instant (frame) of the MD trajectory, 
the atom-atom displacement can be written as AL(r) = 
ALg^^^ + AL^jy^(r). Thermal effects, at least in an approximate 
way, might be added to the simulation of the powder pattern 



[e.g. by introducing a Debye-Waller factor and temperature 
diffuse scattering (Warren, 1990; Beyerlein et al, 2012)]. 
However, for the purpose of studying the static component it 
is more convenient to introduce the concept of time-averaged 
atomic coordinates (TACs): the atomic coordinates of all 
atoms in the studied system are averaged over a suitably large 
number of time frames of the MD trajectory, so that 
AL(r) = ALg^^^. In this way we can get rid of the dynamical 
component of atomic displacement and refer the diffraction 
LPA to the static component only. Once the TACs are known, 
the DSE pattern is easily obtained as explained before. The 
D-PDF pattern is then built according to equations (1) and 
(5). 

If we now move to the system after MD equilibration, 
Fig. 7(a) shows the detail for the 200 peak, as obtained by 
adding all {200} reflections in the cluster. While the right inset 
highlights the trend for the (/zOO) directions in different grains. 
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Figure 7 

Powder patterns of the system of Fig. 1 after MD equilibration, (a) Detail of the 200 peak as built from the D-PDF method of equation (5) [whole 
microstructure (thickest line), single direction of each grain (thin lines)]; (b) comparison between DSE (circles) and D-PDF patterns (line), with 
difference plot below. 
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the plot in the left inset shows an interesting feature of the 
position (^b) of the hOO peaks in the cluster, which changes as 
a function of the size of the corresponding grain, as an effect of 
the strain dependence on the grain size (cf. Fig. 3). 

Fig. 7(b) shows a comparison between patterns for the 
equilibrated system, as obtained by DSE and D-PDF. Small 
but visible differences are expected and suggest interesting 
features of the strain field. While the DSE includes all regions 
of the cluster independently of their order (whether they are 
crystalline or amorphous), the D-PDF refers by definition to 
an underlying crystalline lattice: local atomic displacement is 
allowed (D-PDF position, width and shape) but is in any case 
bounded by the average crystalline framework. The residual in 
Fig. 7(b) is therefore related to a highly strained region, 
somewhere in between a highly distorted crystalline lattice 
and an amorphous phase, which is very likely the grain 



boundary area. Moreover, while the DSE is sensitive to all 
possible correlations within the cluster, i.e. between atoms 
inside each grain as well as between atoms of different grains 
(Leonardi et al., 2013Z?), the D-PDF considers just those 
between atoms inside grains. 

To test this hypothesis, the DSE and the D-PDF patterns 
were calculated again on the same equilibrated cluster but 
after removing, respectively, one, two or three external atomic 
layers from all grains. In this way the grain boundary region is 
progressively eliminated (Fig. 8). 

As shown in Fig. 8(c) for a restricted angular region at 
relatively high q, the match between the DSE and the D-PDF 
patterns is markedly improved by removing surface layers. As 
soon as the first layer is removed, just fully coordinated atoms 
are considered and most of the discrepancy disappears 
(Fig. 8b). Besides zeroing the peaks in the low-^ region of the 





q = 47csme/X, (A ) 

Figure 8 

(a) MD-equilibrated cluster of Fig. 1 after removal of one atomic layer from each grain, (b) DSE pattern for the equilibrated cluster (open circles) and 
after removal of one (upward triangles), two (downward triangles) and three (squares) layers. The line refers to the corresponding patterns calculated 
from the D-PDF approach of equation (5). The difference between the DSE and D-PDF patterns is shown below, (c) A detail for the 531 and 442/600 
peaks. 
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Result of profile fitting of the pattern given by the Debye scattering equation for the studied metallic cluster: DSE data (circles), best fit with pseudo- 
Voigt functions (line) and difference between the two (residual, line below). Corresponding Miller indices are shown in the log-scale plot in the inset. 
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residual, this step eliminates nearly completely the diffuse 
scattering at higher q values. This last detail is visible by 
comparing the inset of Fig. 7(b) and Fig. S(b). 

A progressive removal of layers can thus help the investi- 
gation of the degree of disorder in the grain boundary region, 
the contribution of this region to the strain field in the grains, 
and the corresponding line profile broadening and diffuse 
scattering. The results suggest that the grain boundary region, 
as obtained from the adopted MD equilibration procedure, is 
highly distorted: it is probably not truly crystalline, but it is not 
even a completely amorphous phase (Leonardi et ai, 2013^?). 
Effects on line profiles are clearly visible, even more so the 
smaller the crystalline domains. 



5. D-PDF and r.m.s. strain 

It is interesting to consider again the DSE pattern as an 
'experimental' pattern from a nano-polycrystalline system, 
and make a traditional line profile analysis. As a first step, 
which is ordinary practice in real cases of study, profile fitting 
is used to separate contributions from the different, strongly 
overlapping peaks (Dong & Scardi, 2000). The fitting results 
shown in Fig. 9 are reasonably good, although less satisfactory 
for the peaks at lower q, which are more affected by the grain 
boundary and the grain-grain correlations discussed in the 
previous section. From this analysis (Dong & Scardi, 2000) it is 
straightforward to obtain the plot of Fig. 10, i.e. the logarithm 
of the Fourier transform of the line profiles as a function of q^ 
[actually (q/lnf' for historical reasons] for different pair 
distances L. 



According to Warren and Averbach, the observed data can 
be described as (Warren & Averbach, 1950, 1952; Warren, 
1955, 1959, 1990) 

ln[y,,,(L)] ^ ln[AL(LX,(L, ^)] 

= ln[AL(L)] - 27r'LH6l^(L))(q/27rf. (6) 

Information on the domain size and variance of the strain 
distribution {s\j^i{L)) [or of the displacement distribution, 
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Figure 10 

Warren- Averbach plot: logarithm of the Fourier transform of the peak 
profiles (from profile fitting in Fig. 9) as a function of the square of the 
scattering vector (q/ljT), for a selection of pair distances L (5, 10, 20, 50, 
100 A). Points with the same {qlln) correspond to a given set of Miller 
indices, shown on top for the [hhh] and {/zOO} families (333 and 600 are not 
considered as they overlap with other reflections). 
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{(AL)lj^i) = {6lj^i(L))] can be obtained, respectively, from 
the intercept and slope of the trends in Fig. 10 considered for 
different pair distance values. To properly account for a 
possible dependence on the crystallographic direction, the 
procedure is performed separately for peak profiles belonging 
to different [hkl] families (Warren, 1959), although an analysis 
involving all observed peaks can also be informative. 

Fig. 11 shows the standard deviation ((AL)^^^) as a 
function of L for 111/222/444 and 200/400/800 (333 and 600 
were excluded as they overlap, respectively, with 511 and 442). 
As shown in the original paper by Warren & Averbach (1950), 
all trends should start from the origin: while this happens for 
the [hOO] family, it is not the case for the [hhh] family, which 
gives a trend crossing the abscissa just above L = 20 A. This is 
an artefact at least partly caused by systematic errors in profile 
fitting (e.g. of the 111 peak, see Fig. 9) and by the severe 
overlapping of the peak profiles (especially relevant for the 
weak 444 peak), which are typical problems also in real cases 
of study. However, as expected for the elastic anisotropy of 
copper, the r.m.s. displacement is higher for [hOO] than for 
[hhh], while results from other [hkl] families fall in between 
these two limits. 

Also shown in Fig. 11 is the result from all observed peak 
profiles: this procedure disregards the elastic anisotropy but 
helps in averaging the effects of a non-perfect fitting of the 
peak profiles. 

Several studies propose an interpretation of the trends in 
Fig. 11. According to Adler & Houska (1979), data of Fig. 11 
should obey a simple power law, ((AL)^^^)^/^ = kU^^, with 
-1 < r < 0. The value found for [hOO] is r = -0.44 (2), 
whereas the fit to the data from all peak profiles gives r = 
—0.47 (1). A value around r = —0.5 is considered as typical of 
cold-worked metals, and would be a result of the non-uniform 
strain field of dislocations (Adler & Houska, 1979). 

Before considering the correct trend of {(AL)^^^) from 
the D-PDF analysis, it is worth analysing the results on the 
domain size obtained from the intercepts in Fig. 10. For 



simplicity we consider the [hOO] family, which gave a plausible 
trend in Fig. 11. It can be shown that, in a plot of the size 
coefficients A^^^(L) as a function of L, the intercept with the 
abscissa of the tangent to the curve at L = 0 is a surface- 
weighted mean domain size, (D)^; the second derivative of 
AIj^i(L), instead, is proportional to the length distribution of 
the domain along the scattering vector direction (the so-called 
column length distribution), from which a volume- weighted 
mean domain size {D)y can be calculated (Bertaut, 1950, 1952; 
Warren, 1990). Fig. 12 shows the column length distribution 
along (/zOO); the mean sizes are (D)^ = 48 A and {D)y = 62 A. 

These results can be compared with the values provided by 
the D-PDF analysis, which can be considered 'exact' in that 
they are directly obtained from the known parameters of the 
cluster. The CVF is obtained from the area below the D-PDF 
curves (e.g. those in Fig. 4), and after normalization to the 
grain volume it provides the size coefficients for the given 
grain and [hkl] direction. A^^^(L) is then calculated by aver- 
aging the coefficients over all grains in the cluster, and the 
second derivative provides the column length distribution. 
This procedure, applied to the (hOO) directions, gives (D)^ = 
40 A and {D)y = 59 A, and the column length distribution 
shown in Fig. 12. Despite the quite different procedures 
involved, the column length distributions are remarkably 
similar. The profile fitting/Warren-Averbach method tends to 
overestimate the size, the discrepancy being larger for (D)^ 
than for (Z))y, as the former is more influenced by the shorter 
lengths; these in turn depend more directly on the peak taUs, 
which are less accurately described by the profile fitting. 

Standard deviations of the strain distribution are easily 
calculated numerically for each D-PDF curve, and averaging 
this result over equivalent directions of all grains provides 
{(/^L)lj^i)^'^ for the cluster. The trends for the {/zOO) and {hhh) 
directions are shown in Fig. 13(a), together with the corre- 
sponding results from the Warren-Averbach analysis. 

While the agreement between results on domain size is 
good, the standard deviations of the displacement distribution 
look completely different. Reasons for this discrepancy can be 
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Warren plot: standard deviation of the distribution as a function of the 
atomic pair distance along (/zOO) (squares) and {hhh) (circles), and the 
average over all {hkl) (filled diamonds); the best fit (dash-dot line) refers 
to a power law (see text for details). 




-| I I I I I I I I I I I I I I I I I I I I I |— 
0 10 20 30 40 50 60 70 80 90 100 110 120 130 140 150 

Column length (A) 

Figure 12 

Column length distributions along (/zOO), as obtained from the Warren- 
Averbach analysis of profile data from the DSE pattern (circles) and 
directly from the D-PDF analysis (squares). 
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found in the fine features of the D-PDF curves, and in the 
simphfying assumptions underlying the Warren-Averbach 
method. As shown in Fig. 14, the assumption of a Gaussian 
and symmetrical pj^(AL^j^i) is far from being correct; more- 
over, the shape of the strain distribution changes for different 
pair distances (Fig. 14b). 

However, the non-Gaussian (and asymmetrical) nature of 
p^(AL^i^i) alone cannot explain the large discrepancy in the 
results of Fig. 13. An interesting feature, visible in the log-scale 
plot of Fig. 14(a), is the constant 'background', i.e. the fact that 
the distributions do not fall to zero with increasing distance 
from the expected (perfect crystal) value of atomic pair 
distance L. This constant component is due to the highly 
disordered grain boundary region, which contributes to the 
peak profiles with an atomic displacement effect independent 
of L. This is similar to the effect explained by Warren and 
Averbach for the thermal vibrations, with the important 



difference that the strain component involved here is static, as 
the dynamic component was removed before producing the 
DSE pattern. 

To test this hypothesis, starting from the 'crystallographic' 
system (i.e. before any energy minimization and thermaliza- 
tion procedure), we added increasing levels of random static 
disorder (as in a 'frozen' thermal effect). The D-PDF analysis 
gives the {(AL)^^^) trends in Fig. 15(a), showing a constant 
increase with the equivalent temperature, up to ((AL)^^^) 2:^ 
0.23 for an effect equivalent to the melting point of copper (ca 
1358 K). As shown by the dashed line in Fig. 13, this value 
agrees quite well with the sharp step in the ((AL)^^^) trend. 

A further contribution to the observed strain is of course 
given by the non-uniform strain across grains, and from grain 
to grain, shown in Figs. 2 and 3. These are responsible for the 
L-dependent component of the trend in Fig. 13, and could be 
schematically labelled as strains of II (intergranular) and III 
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Figure 13 

Warren plot comparing results of Fig. 11, obtained from the Warren-Averbach analysis of profile data from the DSE pattern (filled symbols), with the 
standard deviation of the strain distribution directly calculated from the D-PDFs. Results refer to the microstructure of Fig. 1 after MD equilibration {a) 
and after removal of one atomic layer from the surface of all grains in the cluster {b). The full line in {a) indicates the equivalent thermal strain at melting 
temperature. 




Figure 14 

{a) Example of a pair distance distribution (average over the whole cluster for the (/^OO) directions) with the best fits of a Gaussian and of a Lorentzian 
function; {b) corresponding strain distributions for different L values. 
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Figure 15 

Standard deviation of a strain distribution simulated by adding increasing 
levels of random static disorder, with indication of the corresponding 
equivalent temperature. The type II strain effect is also shown (see text 
for details). 

(intragranular) type (Van Houtte, 1993; Hutchings etal, 2005). 
For a simplified estimate of this complex effect, we added a 
pure type II strain to the crystallographic cluster. This was 
done by changing the mean unit-cell parameter of each grain 
according to the values of Fig. 3, resulting in a nearly linear 
increase of the standard deviation with the pair distance 
(Fig. 15). A combination of the two effects, type II strain and 
random static disorder equivalent to 1358 K, explains at least 
qualitatively the trend observed in Fig. 13. The differences are 
due to the simplification of considering a type II strain only, 
whereas the strain also changes inside each grain with a type 
III component. Elastic anisotropy should be taken into 
account too. 

Indeed, as shown in Fig. 13, the standard deviation of the 
atomic displacement distribution for the two extreme direc- 
tions (/zOO) and {hhh) depends on the crystallographic direc- 
tion. Further evidence is provided by Fig. 3(b), where a 
directional dependence of the trend of the unit-cell parameter 



versus the largest thickness is shown. This anisotropic beha- 
viour can be investigated following the idea originally 
proposed by Stokes and Wilson (Sibson, 1980). In Fig. 16 the 
variance ((AL)^^^) for different L values is plotted as a func- 
tion of the invariant for the (cubic) Laue group of copper, 
H = (h^k^ + k^P + Ph^)/(h^ + /c^ + Pf. Excluding the very 
low L values (up to about 20 A), the trend is reasonably linear. 
It should be considered that the first step in the D-PDF is of 
different length for the different hkls (high indices can have 
rather long distances between neighbours). This can explain 
the scatter observed for low L values, where the highly 
disordered grain boundary contribution can play a major role. 

It is also possible to plot the standard deviation {(AL)^^^) 
as a function of H (Fig. 16). The fit is again acceptably hnear, 
but the best regression is obtained for the plot of ((AL)^^/) 
versus H. We could further speculate on this result, to find the 
best functional dependence of ((AL)^^^) on //, but the main 
conclusion is that this correlation is a consequence of the 
elastic anisotropy of the metal. A similar effect is also 
observed when dislocations are responsible for line broad- 
ening and can be expressed via the contrast factor concept 
(Wilkens, 1969, 1970; Martinez- Garcia et aL, 2009), but in our 
case, quite evidently, no dislocations are present. This suggests 
great caution is required in univocally attributing line profile 
broadening anisotropy to dislocations, which, quite evidently, 
are just one possible source for this directional effect. 

The analysis presented in this section can be repeated on 
the DSE patterns of Fig. 8. After removal of one surface layer 
from each grain (corresponding to removal of all atoms with 
coordination lower than 12), all D-PDFs are in general much 
narrower than those of the cluster completely equilibrated by 
MD. In addition (see Fig. 17), the removal of a surface layer 
drastically reduces the extension of the distribution tails and 
of the constant component, which was interpreted as a 'frozen' 
thermal-like effect of the disordered grain boundary region. 
An additional layer removal further extends this effect. 

The standard deviation of the atomic displacement is 
consequently smaller than in the starting cluster. The observed 
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Figure 16 

Variance of the strain distribution as a function of the invariant form H, for different values of the pair distance, L. Results refer to the MD-equilibrated 
system before (a) and after (b) removal of one atomic layer from all grains. 
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Figure 17 

Example of a distribution from Fig. 14(a) for the MD-equilibrated system 
(filled circles) and after removal of one atomic layer from all grains (open 
circles). 

trend in the Warren plot (Fig. 13b) is much more 'regular', i.e. 
as expected from the Warren-Averbach model (Warren & 
Averbach, 1950), as a consequence of the much reduced 
(nearly eliminated) constant contribution of the grain 
boundary. This further confirms the role of the grain boundary 
as a region with high disorder but which still contributes to the 
coherent scattering of the crystalline grains: not all scattering 
from the grain boundary is therefore diffuse scattering. As a 
further support to this interpretation. Fig. 13(b) shows that the 
WA analysis on the DSE pattern after removal of one surface 
layer from each grain leads to results much closer to those 
from the D-PDF analysis than for the starting cluster (cf. 
Fig. 13a). Finally, it is also worth noting that after removal of 
one layer the anisotropy of the line broadening is more 
evident, even if the total strain is lower than in the original 
cluster. As shown in Fig. 16(b), the plot of the variance as a 
function of H gives linear correlations for any L value with a 
lower data scatter than in Fig. 16(a). 

6. Conclusions 

The concept of the directional pair distribution function has 
been introduced to support a better understanding of line 
broadening effects in the diffraction peak profiles from nano- 
polycrystalline microstructures. The new concept, illustrated 
for a simple system made up of (nearly equiaxed) nanocrys- 
taUine Cu grains equilibrated by MD, can easily be applied to 
any simulated microstructure, including lattice defects and 
crystalline domains of any shape. 

The D-PDF analysis shows in detail how the local atomic 
displacement influences line profiles, also taking into account 
the anisotropy of the strain distribution. Most importantly, the 
D-PDF approach leads to results equivalent to a traditional 
line profile analysis based on a Fourier formalism, such as the 
method of Warren and Averbach, thus providing a direct 
possibility to understand the meaning of LPA results in terms 
of an atomistic model of the microstructure. This capability 
can be especially useful for a correct interpretation of 



experimental LPA results and can provide further insights into 
the analysis of diffraction phenomena. 
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